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Abstract 

We consider a one-dimensional totally asymmetric nearest-neighbor 
zero-range process with site-dependent jump-rates — an environment. 
For each environment p we prove that the set of all invariant mea- 
sures is the convex hull of a set of product measures with geometric 
marginals. As a consequence we show that for environments p sat- 
isfying certain asymptotic property, there are no invariant measures 
concentrating on configurations with density bigger than p* (p) , a crit- 
ical value. If p*{p) is finite we say that there is phase-transition on 
the density. In this case we prove that if the initial configuration has 
asymptotic density strictly above p* (p) , then the process converges to 
the maximal invariant measure. 

AMS 1991 subject classifications. Primary 60K35; Secondary 82C22. 
Key words and Phrases. Zero-range; Random rates; invariant mea- 
sures; Convergence to the maximal invariant measure 



1 Introduction 

The interest on the behavior of interacting particle systems in ran- 
dom environment has grown recently: Benjamini, Ferrari and Landim 
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(1996), Evans (1996) and Ferrari and Krug (1996) observed the ex- 
istence of phase transition in these models; Benjamini, Ferrari and 
Landim (1996), Krug and Seppalainen (1999) and Koukkous (1999) 
investigated the hydrodynamic behavior of conservative processes in 
random environments; Landim (1996) and Bahadoran (1998) consid- 
ered the same problem for non-homogeneous asymmetric attractive 
processes; Gielis, Koukkous and Landim (1998) deduced the equi- 
librium fluctuations of a symmetric zero range process in a random 
environment. 

In this article we consider a one-dimensional, totally asymmet- 
ric, nearest-neighbor zero-range process in a non-homogeneous envi- 
ronment. The evolution can be informally described as follows. Fix 
c G (0, 1) and provide each site x of Z with a rate function px £ [c, 1]. 
If there is at least one particle at some site x, one of these particles 
jumps to a; + 1 at rate px- A rate configuration p = {px : x € Z) is 
called an environment and a measure m on the set of possible envi- 
ronments a random environment. 

Benjamini, Ferrari and Landim (1996) and Evans (1996) for an 
asymmetric exclusion process with rates associated to the particles 
— which is isomorphic to a zero range process with rates associated 
to the sites — and Ferrari and Krug (1996) for the model considered 
here, proved the existence of a phase transition in the density. More 
precisely, they proved that, under certain conditions on the distribu- 



tion m, specified in Theorem |2.4| , there exists a finite critical value 



p* such that for m-almost-all p there are no product invariant mea- 
sures for the process with rates p concentrating on configurations with 
asymptotic density bigger than p* and that there are product invari- 
ant measures concentrating on configurations with asymptotic density 
smaller than or equal to p* . (The density of a configuration is essen- 
tially the average number of particles per site and is defined in ^ 
below) . 

Our first result is that the set of extremal invariant measures for 
the process with fixed environment p = {px : x G Z) is the set {vp^v '■ 
V < Px,yx}, where is the product measure on with marginals 

^P,.{e:^(x) = A:}=(f)Vl-f). (1) 

^ Px ' ^ Px ' 

The above result does not surprise specialists in queuing theory. In 
fact we are dealing with an infinite series of M/M/1 queues with service 
rate px at queue x. The value v can be interpreted as the arrival rate 
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at "queue" — oo. Since Burke's theorem (see Kelly (1979) or Theorem 
7.1 in Ferrari (1992) for instance) guarantees that in equilibrium the 
departure process of a M/M/1 queue is the same as the arrival process 
(both Poisson of rate v), there is an invariant measure for each arrival 
rate v strictly smaller than all service rates. 

Assume c = inf^Px and that the following limits exist. For v < c, 

1 ° /■ 1 ° V 
Rip,v) := \im - y~] / fp,„((i^) ^(x) = lim - . 

x=—n+l x=—n+l 

(2) 

We interpret R{p, v) as the global expected left density per site of 
the configurations distributed according to I'p^y. A consequence of the 
existence of the limits, as we will explain later, is that for all v < c, Vp^^ 
concentrates in configurations with asymptotic left density R{p,v): 

up^ J lim- V C{x) = R{p,v)) = 1 . (3) 

Vn— >oo n ^ — ' / 

x=—n+l 

It is easy to prove that R{p, v) is a strictly convex increasing function 
of V, hence the limit 

p*{p) := liniR{p,v) (4) 

is well defined (but may be infinite). In the sequel we assume p*(p) < 
oo. We do not assume the existence of the limit in (^) for v = c, nor 
the fp^c almost sure convergence of the density. 

Our second and main result states that under the condition p* (p) < 
oo, initial measures concentrating on configurations with asymptotic 
left density strictly bigger than p*{p) converge towards the maximal 
invariant measure Vp^c- We do not know in general if this measure 
concentrates on configurations with density. But if the limit R{p, c) 
of (HI) exists, equals p*{p) and is finite our result says that the process 
starts with global density strictly above p* (p) and converges to a mea- 
sure with density p*{p)- This behavior is remarkable as the process is 
conservative, i.e. the total number of particles is conserved, but in the 
above limit "looses mass" . Informally speaking, what happens is that 
many clients remain trapped in far away slow servers. More precisely, 
denoting by Sp{t) the semigroup of the process, we first show that for 
any initial measure z^, all weak limits of the sequence {i'Sp{t), t > 0} 
are dominated, in the natural partial order, by Vp^c- We then show 
that if is a measure concentrated on configurations with asymptotic 
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left density strictly greater than p*{p), all weak limits of vSp{t) dom- 
inate Vp^c- Surprisingly enough, the proof of the second statement is 
much more demanding than the proof of the first one. 

It follows from the two previous results that the domain of at- 
traction of c includes all measures with asymptotic density strictly 
above p*{p)- It remains an open question to describe the domain of at- 
traction of a product invariant measure h'p^y for <^ v <^ c or to show 
the convergence to Up^c of initial measures with asymptotic density 
p*{p). 

Our results hold m-a.s. for measures m concentrating on environ- 
ments satisfying (^). 

The paper is organized as follows: in Section 2 we introduce the 
notation and state the main results. In Section 3 we characterize 
the set of invariant measures and show that the maximal invariant 
measure dominates all the weak limits of the process. In Section 4 
we obtain the asymptotic velocity of a second class particle for the 
zero-range process in a non homogeneous environment and use this 
result to prove the main theorem. 

Many of our results are based on standard coupling arguments. We 
assume the reader familiar with this technique described in Section 1 
of Chapter 2 of Liggett (1985). 

2 Notation and Results 

Fix < c < 1 and consider a sequence {px)x<az taking values in [c, 1] 
such that c = mf^Px- We consider a totally asymmetric zero-range 
process in the environment p. This is a Markov process that can be 
informally described as follows. We initially distribute particles on 
the lattice Z . If there is at least one particle at some site x, then 
at rate px one of them jumps to site x + 1. To construct a Markov 
process r/j on X = corresponding to the above description, let 
^x{t) (x G Z) be a collection of independent Poisson processes such 
that for all a; G Z, E{Nx{t)) = Pxt. The evolution of r]t is now given 
by the following rule: if the Poisson process Nx{-) jumps at time t and 
r]t- (x) > then one particle is moved from x to x + 1 at that time. To 
see that the process is well defined by this rule, just note that in any 
time interval [0, t] for any x there exists with probability 1 a y < x 
such that Ny{t) = 0. Hence the value of r]t{x) depends only on the 
initial configuration and on a finite number of jumps. 
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The generator Lp of this process, defined by Lpf{r]) = d'E[f{r]t) \ % 
v]/dt 



, acts on cylinder functions / as follows: 

t=o 



In the above formula r]^ = r] — dx + Ox+i; where dy stands for a con- 
figuration with just one particle at y and addition of configurations is 
performed componentwise. 

We denote by {Sp{t), t > 0} the semigroup associated to the gen- 
erator Lp, i.e. Sp{t)f{r]) = E[/(7/j) \r]o = rj] and by Ip the set of 
invariant measures of r]t (the Markov process with generator Lp). Let 
t; be a real number such that < v < px for all x. Then a standard 
calculation (first observed by Jackson (1957) for the finite case) shows 
that the product measure Vp^y with marginals given by ([l|) is an invari- 
ant measure for the process. Benjamini, Ferrari and Landim (1996) 
raised the question of whether or not there exist invariant measures 
which are not convex combinations of the fp,^'s. In Section 3 we prove 
the following theorem which, combined with Theorem 12.2 in Dynkin 
(1978) (which states that the set of extremal invariant measures is 
the convex hull of the set of invariant measures) , gives a negative an- 
swer to that question. In its statement we denote by {Ip)e the set of 
extremal invariant measures for the process. 

Theorem 2.1 Let p be an arbitrary environment then 
{Ip)e = {I'p^v ■ V <px, yx e Z}. 

In this theorem the range of the parameter v may be either [0, c) or 
[0, c] - when Px = c for some x or px > c for all x, respectively. In the 
first case to prove the theorem we only need to follow the proof of The- 
orem 1.11 in Andjel (1982), but in the second complementary 
argument is needed. In both cases the proof relies on the standard 
partial order for probability measures on X. To define it, first say 
that r/ < ^ if r]{x) < ^(x) for all x G Z. Then say that a real valued 
function / defined on X is increasing if 7/ < ^ implies that f{r]) < /(^). 
Finally if /i and are two probability measures on X, say that n <v 
if J fd^ < J fdv for all bounded increasing cylinder functions /. In 
this case we say that v dominates /i. The complementary argument 
alluded above depends on the following proposition: 
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Proposition 2.2 Assume that p is an environment such that 



Px > c for all X G Z and liminf = c, (6) 



and let v he an arbitrary probability measure on X. Then the set of 
measures {i'Sp{t) : t > 0} is tight and its weak limits as t goes to 
infinity are bounded above by Vp^c- 



An immediate corollary of Proposition 2.2 is that under (y) all 
invariant measures are dominated by Vp^c- 

To state our main result let r/ be an element of X and consider 

1 ° 

^(77) = liminf— > 77(2;) , 

- 1 ° 

D(r]) = limsup— ; 



n— >oo , , 



the lower, respectively upper asymptotic left density of r]. If both limits 
are equal to a we say that r] has left density a and write D{ri) = a. 

Assume that p is an environment for which the limits defined in (|2|) 
exist. Then, by Kolmogorov's law of large numbers (see e.g. Shiryayev 
(1984), Theorem 2 p. 364) fp ^, concentrates on configurations with 
left density R{p,v): 

Up,,{rjeX:Di7]) = R{p,v)} = l (7) 

for all V < c. 

The values assumed by R{p, v) for v < c are crucial for the char- 
acterization of the set of invariant measures for the process with rates 
p. If liuiy^c R{p, v) = 00, then the range of allowed densities is [0, 00) 
or [0, cxd]. The first case occurs when p^ = c for some x. In this case 
t'p.t, is defined for any v < c, but not for v = c. Moreover, since R{p, •) 
is continuous and increases to 00 as u — > c, then for all p € [0, 00) 
there exists v = v{p,p) such that Vp^vlv G X : D{rj) = p} = 1. The 
second case occurs when px > c for all x. In this case fp^c is well 
defined and concentrates on configurations with infinite asymptotic 
left density, and for any p G [0, 00] there exists v = v{p, p) such that 
Vp.,{r] £ X : D{7]) = p} = I. 

If limy^cRip,v) = p*{p) < 00 and px > c for all x, the measure 



Up^c is well defined and Theorem 2A tells us that there are no invariant 
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measures bigger than fp^c- Our next theorem describes what happens 
in this case when one starts with a density strictly bigger than p*{p). 
This is our main result. 

Theorem 2.3 Letp he an environment satisfying ^ such that p*{p) < 
oo and rj be a configuration such that Div) > P*{p)- Then 

lim S„Sp{t) = fp^c, 

t — >oo 

where 5n is the measure giving weight one to the configuration rj. 



As a corollary to Theorem 2.2 we obtain the asymptotic behavior 



of the system when the environment is randomly chosen. Let m be 
the distribution of a stationary ergodic sequence p on [c, 1] such that 
'm{{p : pq = c}) = 0, m{{p : c < pq < c + e}) > for all e > 0. 
The f = J m{dp) f Vp^v{dri) f [ri) is an 

ergodic distribution on X and, by the Ergodic Theorem, for all w < c 
and for m-almost all p, the asymptotic density exists Vp^y a.s. and is 
equal to: 

[ V 

R(v) = / m(dp). 

J Po-v 

Let p* := limt,_»c -R(i') and assume p* < oo. In this case for m- 
almost all environment p any invariant measure for Lp is dominated 
by fp,c- The following theorem concerns the behavior of the process 
when the initial measure concentrates on configurations with density 
strictly higher than p* . 

Theorem 2.4 Let m be the distribution of a stationary ergodic se- 
quence p = {px)xez on (c, 1] such that m{{p : c < po < c + e}) > for 
all e > and for which p* < oo. Let v he a measure for which v a.s. 
Diji) 'i'S strictly bigger than p* . Then, for m-almost all p 

lim uSpit) = Up^c- 



3 Domination and Invariant measures 



In this section we prove Proposition p.2| and Theorem ^ . 

Proof of Proposition 2.2. Fix an arbitrary site y and let x„ be a 

decreasing sequence such that xi < y, px„ < pz ior Xn < z < y and 
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Pxn decreases to c. The existence of such a sequence is guaranteed by 
(|6|). Consider a process on pj^^'""^^'"''^^ with generator given by: 

Lp,nf{ri) = l{^(^)>0}p.[/(r/^)-/(r/)] 

+ Px.J/(r? + fx„+i)-/(r/)] (8) 
+ l{r?(y) >0}p,[/(r/-i),)-/(77)] . 

Let Sp^n be the semigroup associated to this process and for an ar- 
bitrary probabihty measure v let m„ be its projection on 
Standard couphng arguments show that 

The couphng of the two processes is done using the same Poisson 
processes Nx{t) defined in Section 2. The reason why the domination 
holds is that for the process Sp^niJ^), each time the Poisson process 
Nx„ (t) jumps, a new particle appears in while the same happens 

for the process Sp{t) only when there is at least a particle in the site 

The process with generator Lp^„ is irreducible and has a countable 
state space, moreover a simple computation shows that the product 
measure fin^p with marginals given by 

^^^,p{^:^{z)=k}=(l-P^)(P^y, 

where x„ < z < y, is invariant for the process. Therefore VnSn,p{t) 
converges to ^n,p and any weak limit point of {vSp{t))n is bounded 
above by /in,p- Since as n goes to infinity the marginals of ^n.p converge 
to the marginals of fc,p the proposition is proved. ■ 



Proof of Theorem 2.1. Since only the final step of the proof is 
different from the proof of Theorem 1.11 in Andjel (1982) (in which the 
set of all invariant measures is characterized for a family of asymmetric 
zero-range process) we refer the reader to that paper. Exactly as there 
one proves that if Vp is an extremal invariant measure then for each 
V < inlxPx either fp < Vy^p or > v^^p. This implies that either 
i^p = T^v,p for some v ot Vp > v^^p for all v. The latter case cannot 
occur if there exists x such that px = \niyPy because this would imply 
that Vp{r] : r]{x) > fc} = 1 for all k. Therefore Vp > Vc,p and either 
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liminf^^_oo Px = infpy or liminf2.^oo Px = vnipy In the first of 
these cases, Proposition allows us to conclude immediately. In 
the second case we argue by contradiction: let ? be a probability 
measure on x admitting as first marginal and second marginal 
Vp and fc^p respectively and such that ^'{(??, ■ > ^} = \. Consider 
the standard coupled process with initial measure V. Denote by S{t) 
the semigroup associated to this process and assume that for some x, 
^{{ViO '■ > Ci^)} > 0- Suppose it exists k and Z in N \ {0} such 
that i/{r/(x) = k + l,(,{x) = A;} = ei > then at any time 6 > one 
can find a £2 > such that 

uS{6){ri{x) = k+l-l,^{x) = k-l} = u{r]{x) = k+l-l,C{x) = k-1} = £3. 

To see that, one has just to control the arrivals and departures of 
particles on sites x — 1 and x which are given by exponential clocks. 
By induction it follows that for all t > 

i?5(t){(7?,e) >e(x) =0} >0. 

Hence 

Vp {r] : r]{x) > 0} > i^cp {r] ■ r]{x) > 0} = — . 

Px 

Pick y > x and such that py < Pxi^pii] ■ r]{x) > 0}. Then let = 
^^z=x+i'ni^)- Now a simple calculation shows that / Lpf{ri)di'p{rj) > 
contradicting the invariance of fp. ■ 



Remark: Proofs of Theorem 2.1 and Proposition 2.2 can easily be 



extended to a larger class of one-dimensional nearest-neighbors asym- 
metric zero range processes in non-homogeneous environment. In 
these systems a particle at site x on configuration tj jumps at rate 
Pxdi^ix)) to site X + 1, where : N — > [0,oo) is a non-decreasing 
bounded function such that g{0) = 0. 



4 Convergence 



We prove in this section Theorem 2.2. Fix a measure u on X con- 
centrated on configurations with lower asymptotic left density strictly 
greater than p*{p)- Let z> be a weak limit of vSp{t). Proposition 2.2 



shows that v is dominated by Vp^c- Lemma [4.1| below implies that 
D dominates Vp^^ for all v < c. This finishes the proof of Theorem 



2.3 because : < u < c} is an increasing sequence of measures 
converging to Vp^c- 
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Denote {Sp{t) : t > 0} the semigroup corresponding to the cou- 
phng between two versions of the process with (possibly) different 
initial configurations, by using the same Poisson processes {Nx{t) : 
X G Z) in its construction. 

Lemma 4.1 Let p be an environment satisfying ^ and such that 
p*{p) < oo and ( o- configuration with lower asymptotic left density 
D_{C,) > p*{p). Then for any v < c, 

lim{S^ X Up,,)Sp{t)[{7],0 : r/(x) < ^(x)} = (9) 

for all X in 7j. 

The proof of this lemma requires the following result. It states that 
for each v < c the asymptotic velocity of a second class particle in the 
zero-range process in the environment p under the invariant measure 
fp^„ is strictly positive. 

Fix a starting site z and consider a coupled zero range process 
with initial condition (r],rj + dz) and semigroup Sp{t). Under the 
coupled dynamics the number of sites where the marginals differ does 
not increase in time. Let be the site where the marginals differ 
at time t. We can think that X^ stands for the position of a "second 
class particle" . Indeed, if the second class particle is at x at time t it 
jumps to x + 1 at rate px l{r]t{x) = 0}. In other words, the second 
class particle jumps only if there is no other particle at the site where 
it is. 

For an environment p and a probability measure on X, denote 
by Pjy the measure on D(R_|_, X) induced by i/ and the Markov process 
with generator Lp defined in (^). In the next lemma we write P(i/,z) 
for a coupled process whose initial configuration is (r] , rj + D z) , with 
r] distributed according to v. Since R{p, v) is convex and strictly 
increasing 

^{p,v):=[R\p,v)\-^ (10) 

exists in a dense subset of (0,c). In the sequel we abuse notation by 
not writing integer parts where necessary. 

Lemma 4.2 Letp be an environment for which the limits in exist 
for V < c. Pick v € (0, c) such that ^{p, v) exists. Then, 



—at 



^limP(,^„,_,j)(|^-(7(p,t;)-a)| >e) =0, (11) 
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for all e > if a > ^{p, v) and 



lij^PK,„,-a4^>0) =1 (12) 



ifa< jip,v). 



Remark. The more complete result when the starting point a is 
greater than 'y{p, v) comes from the fact that in our hypothesis we 
have only the asymptotic left limits (0). If the limits (|2|) hold for both 



sides, then (11) is valid for all a. 

Proof: Note that it suffices to prove ([lT|), since ( p!2[) follows from ( [III) 
because < Xf for all t > if x < y and because (|l^) does not 
depend on the environment to the right of the origin. For u < w < c 
let 

^p,u,w be the product measure on X x X whose first marginal is 
equal to Up^u, whose second marginal is equal to Up^^ and which is 
concentrated above the diagonal: i'p,u,w{{'>]jO '■ ^ < '^j = 1- Denote 
by {rjt,(,t) the coupled Markov process starting from 9p^u,w 

Denote by Ct the difference — rjt and observe that the C-particles 
evolve as second class particles in the sense that a C-particle jumps 
from X to rr + 1 at rate 

[l{7]{x)+ax) > 1}-1{t]{x) > 1}]; 

that is, when there are no r] particles present. In this case we say 
that the ry particles have priority over the C, particles. We label the 
C-particles at time in the following way. Without losing much we 
can assume that there is a C particle at site (integer part of) —at. 
The measure conditioned on this event is absolutely continuous with 
respect to i'p^u,w, an this will be enough for our purposes, as we shall 
only use laws of large numbers. Call particle this particle, and 
complete the labeling in such a way that a particle with label j is at 
the same site or at the left of a particle with label k j < k. Denote 
by the position at time t of the particle labeled j. By construction, 
we have • • • < Yq'^ < Yq = —at <Yq<---. We let the second class 
particles evolve in a way to preserve this order. To keep track of the 
densities involved in the definition we call Yf^'^ = Y^". 

Consider now a single second class particle for the rj process ini- 
tially at the position of —at. This is obtained by considering the cou- 
pled initial condition (ry,ry + i)_at)- Denote the position of the single 
second class particle at time t by X^ (for u = v, this has the same law 
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as the particle denoted by Xf'^ in the statement of the proposition). 
Since Yq'^ = Xq, in the coupled evolution obtained by using the same 
Poisson processes {Nj.{t)) we have Y^'^ < Xf for all t. Indeed, in this 
coupling Y^' for /c > have priority over Y^'^ while those particles 
have no priority over Xf. Similarly, consider a second class particle 
for the ^ process and denote it Xf . Since Y^ for k < Q have priority 
over X]" but not over l^"'*", X]^ < Y^'^ . Hence, for < u < w < c, 

Xf<Yr<X^, (13) 

^i>p,u,-a, almost surely. 

Denote by j/, J}^"^ and J"^ the total number of ry, ^ and C, particles 
that jumped from —at to —at + 1 before time t. In particular, = 
j}^"^ — Jt ■ By Burke's theorem, the number of T/-particles (resp. ^- 
particles) that jump from —at to —at + 1 is a Poisson process of 
parameter u (resp. w). Hence the number of ^-particles that jump 
from —at to — at + 1 in the interval [0, t\ is the difference of two Poisson 
processes and satisfies the law of large numbers: 

t2 t1+2 _ Tl 

lim — = lim = w — u 

t— +00 t t— >oo t 

in Ppp „ ,„ probability. On the other hand, for every t > 0, 

x=—at+\ x=—at-\-\ x=—at+\ 

where \At\ < Ct(Xt^''^)- Note that Q{Yi^''^) is stochastically bounded 
above by a geometric random variable of parameter w/c. Therefore 
\At/t\ converges to in Ppp „ ^, probability as t goes to infinity. As 
in the proof of Theorem 12.1 of Ferrari (1992), it follows from the 
previous equation and the law of large numbers for rjt and — that are 
distributed according to product (invariant) measures with densities 
R{p, u) and R{p, w) respectively — that for u and w strictly smaller 
than c, 

hm-^ h a = ^ 7 (14) 

tloo t R\p, w) — R{p, u) 

in Ppp „ „ probability. Notice that we used here the fact that a > 
7(p, f). In this case Y^'^ /t < and the previous sums refer only to 
negative sites. Hence, from ( p^ we have 

. Xr w-u . X}" 

lim^ + a < — ^ — < lim— ^ + a (15 

tloo t R{p, W) - R{p, U) f\oo t 
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in Ppp „ ^, probability. Fixing w = v and taking the limit n — > u and 
then fixing u = v and taking the limit w ^ v and taking account of 



the differentiability of in we get (11) and (12). 



We are now in a position to prove Lemma ^ 



Proof of lemma 4.1 The proof is performed via coupling. We start 
with two different initial configurations r] and ^ with marginal dis- 
tributions u and fp_„, respectively. Hence r] has lower asymptotic 
density bigger than p*{p) and ^ has asymptotic density R{p,v). We 
use the same Poisson processes for both processes and call {r]t,S,t) the 
coupled process. The configurations t] and ^ are in principle not or- 
dered: there are (possibly an infinite number of) sites z such that 
{rj{z) — ^(z))^ > and (possibly an infinite number of) sites y such 
that (r?(y) — > 0- We say that we have ry^ discrepancies in the 

first case and ^r/ discrepancies in the second. The number of coupled 
particles at site x at time t is given by 

^t{x) ■.= mm{r]t{x),^t{x)} (16) 

The ^ particles move as regular (first class) zero range particles. There 
is at most one type of discrepancy at each site at time zero. Discrep- 
ancies of both types move as second class particles with respect to the 
already coupled particles. When a 77^ discrepancy jumps to a site z 
occupied by at least one ^ry discrepancy, the 77^ discrepancy and one of 
the ^r/ discrepancies at z coalesce into a coupled ^ particle in z. The 
coupled particle behaves from this moment on as a regular (first class) 
particle. The same is true when the roles of ^ and 77 are reversed. 

The above description of the evolution implies in particular that a 
tagged discrepancy can not go through a region occupied by the other 
type of discrepancies. 

We will choose a negative site y such that the jump rate from 
y — 1 to y is close to c. Then we follow the ^r/ discrepancies belonging 
to two disjoint regions of Z at time and give upper bounds on the 
probability of finding them at y at time t. 

Roughly speaking, a ^r/ discrepancy at y cannot come from a region 
"close" to y because we prove that there is a minimum positive velocity 
for the ^7] discrepancies to go. This velocity is given by the velocity of a 
second class particle under Vp^^. On the other hand, the ^tj discrepancy 
cannot come from a region "far" from y because due to the difference 
of densities, a lot of r/^ discrepancies will be between it and y and 
hence they must pass site y — 1 before it. But since we have chosen a 
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small rate for this site, a traffic rush will prevent them to pass. With 
this idea in mind, we have to choose the "close" and "far" regions and 
the value of the rate at y — 1. 

Fix V < c such that R{p, •) is differentiable in v. Let 7 = ^{p,v), 
the (strictly positive) asymptotic speed of a second class particle under 
Vp^y in the sense of (|lO| ) . Denote by f3 the difference between the lower 
asymptotic density of ry and Rip, v) : 



P = P(p,v) = liminf— [r](x) — R(p,v)] . 

n^oo n ^ — ' 

x=—n+l 

For reasons that will become clear later (cf. display (|27|)), we let 

R'{p,v){c-v) 

b = b{p,v) = " < 1, 

p*{p) - R(P,v) 

by the convexity of R; recall that p* [p) = lim^^c 

Rip,v). With this 

choice. 

This allows us to fix e = e{v) satisfying 

< e{v) < Pjb - c + v . 
Finally, choose a negative site y = y{v) such that 

Py-i <c + e. (18) 

We shall prove that 

hm {ly X Up^,)Sp(t)\{v,0 ■■ r]{y) < Ciy)} = . (19) 

t — >oc K. ) 

We can order the ^r/ discrepancies and assume without loss of 
generality that the order is preserved in future times as we did in 
Lemma Of course some of the discrepancies will disappear. Let 

= Zi{^,r]) the positions of the ordered ^r] discrepancies at time 
t with the convention that = 00 if the corresponding discrepancy 
coalesced with a rjS^ one giving place to a ^ coupled particle. Let 

:= {a ^7] discrepancy in the box [y — (t^b), y] at time (20) 
has moved to site y at time t} 

:= Uk{z^G [y-{t^b),y], Z^ = y] 
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where 6 := (l + 6)/2 G (6,1). Hence 



P(A,t(^,6) < P [mm [z^:Z^e [y - itjb),y]] < y) (21) 

We wish to give and upper bound to the event in the right hand 
side above. To do so we consider the coupled (r?,0 process and the ^ 
process to which we add a unique second class particle at y — (t'yb), 
evolving together with jumps occurring at times given by the same 
Poisson processes. We denote by Xf~^"'^ the position of the second 
class particle at time t. If the second class particle has reached y+1 no 
later than time t, then there exists an increasing sequence of random 
times < Ty_(^^^^ < Ty_(^^^^^-^ < ... < Ty such that at each of these 
times the corresponding site has been emptied of its ^ particles. But 
this implies that all the ^rj discrepancies which at time were in the 
interval [y — {tji)) , y] have disappeared or are strictly to the right of 
y. Therefore: 

F{A,AV,0) < P{Xr'^'<y). (22) 



By (|l^) this tends to as t tends to infinity because 76 < 7. The 
above argument is independent of the value of Py^i. 

It now suffices to check that the probability that a ^r] discrepancy, 
to the left of y — (t'jb) at time reaches y no later than time t, tends 
to as t tends to infinity. Let 

:= {a ^r] discrepancy in {—oo,y — (t^b)] at time (23) 

has moved to site y at time t} 
:= Ufc {4 e (-00, y - (t^b)], = y} 

Call W^{ri, C) the positions of the rjS, discrepancies at time t, Wq being 
the first t]^ discrepancy to the left of the origin. As before set = 00 
if the kth. discrepancy coalesced with a £,rj one before t. 

Since a ^ry discrepancy cannot cross over an r]^ discrepancy, 

\z<y—'ytb \^=z 

C {If - 1} > t^(3b} (24) 
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where if and // are the number of rj, respectively particles jumping 
from y — 1 to y in the interval [0, t]. Since 



^o(x)) > t7/36 




Ux)) < t-f(3b } , (25) 



to bound P{B^^tiViO) it suffices to bound the probabilities of the sets 
on the right hand sides of (|2j) and (p5[). For (^) we have 

P(/2 - II > t^ph) < P(iV,^+^ - N^: > t-f(3b), (26) 

where is a Poisson process of parameter a. The above inequality 
holds because the r/-particles jump from y — 1 to y at rate not greater 
than Py-i, which is by construction less than or equal to c + e. On 
the other hand, by Burke's theorem, the number of jumps from y — 1 
to y for the ^-particles is a Poisson process of rate v. By the law of 
large numbers for the Poisson processes, we have 

hm ^(iVf+^ - iV,^) = c-v + e < (3^b, (27) 

t— >oo t 



because we chose e < 7/36 — c + v. Hence (|26D goes to zero as t — oo. 

On the other hand, the probability of the set in the right hand side 
of (H) is 

P( sup _J^(??o(x)-eo(x)) <t7/36) (28) 

By the ergodicity of ^ and the fact that r/ has left density, with prob- 
ability one: 

1 

hm - V (,^o(x)-eo(x)) = 7/36 > 7/36, (29) 

t^co t ^ — ' 

x=y—t'yb 

by the way we chose 6. This implies that ( p8[ ) goes to zero as t — > oo. 
This proves 

To deduce the statement of the lemma from ([T9|) we need the 
following lemma which says that if there exists a subsequence of times 
giving positive probability to a cylinder set, then any other cylinder 
set obtained by moving one particle to the right has the same property. 
These hues follow Andjel (1982). 
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Lemma 4.3 Let f be the following cylinder function on x N^. 

fir], = l{r?(x) = = ^» -.xeA} (30) 

for some finite A C 1^ and arbitrary configurations fi,C ^ N^. Let 
z ^"L be an arbitrary site; define as 

r(r?,0 = l{r/(x) = t{x), i{x) = eix) -.xeA} (31) 

Let {r]t,Ct) be the coupled process starting from an arbitrary measure. 
Then 

limsupE/(r/t,^() > implies limsupE/^(r?(, ^i) > (32) 

t— >oo t— >oo 

Proof. Let A = {x : x + 1 G vl}. Since if z ^ Au A implies that 
= /^(^)0 (and hence for these z the lemma is trivial), we fix 
a z (z AU A. Assume that tn is a sequence of times such that 

hm E/(r?t„,eiJ = c>0 (33) 

n—>oo 

Fix a time s (equal to one, for instance) and consider the event B.„ = 
{N,{tn + s)-N,{tn) = 1, N^{tn + s)-N^itn) = for X € AUA\{z}}. 
That is the event "exactly one Poisson event occurs for z in the interval 
[tm tn + s) and no events occur for the other sites in A U ^ in the same 
time interval". Then 

E/(r?*„,6jP(B„) < E/^(,?i„+„et„+s) (34) 

Since the probability of is independent of n and positive, this 
proves the lemma. ■ 
We continue with the proof of Lemma 4.1. Take an arbitrary y 
satisfying ([l9|). Consider the coupled process starting with the mea- 
sure (z/ X Up^^). By Proposition |2]^ both {r]t{y — l),^t{y — 1)) and 
{'r]t{y),(,t{y)) are tight sequences. Hence there exists a K depending 
on py^i,py such that 

hm P(7?ij2/-1) <Cij2/-l)) >0 (35) 

n^oo 

implies 

hm P(r?j„(y-1) < 6„(y-l), ^tAV'^) < K, mAv) <K)>0 (36) 
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Now we apply Lemma |4.3| : move first the (at most) K rj particles from 
y, then the (at most) K extra ^ particles from y — 1 to y to obtain 
that ( |36[ ) implies 

hm P(r?t,(y) <et;(2/)) >0 (37) 

n^oo 

for some subsequence (i^), in contradiction with (|19|). With the same 
argument we can go to 2; = y — 2, y — 3, . . .. This proves that (|9|) holds 
for all X < ?/ for y satisfying (p^. On the other hand, the marginal 
law of the coupled process at x does not depend on the value of py for 
y > X. Hence, we can assume (|l8|) for y > x + 2 and obtain the result 
for all 2; S Z. This argument works because when we modify Py-i we 
change the process only to the right of y — 1, maintaining the values 
R{p, v) and ^{p, v) unaltered, as they are asymptotic left values. For 
this reason we can use the same e in (p^. ■ 
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